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Abstract. Holder estimates for spatial second derivatives are proved for solutions of 
fully nonlinear parabolic equations in two space variables. Related techniques extend 
the regularity theory for fully nonlinear parabolic equations in higher dimensions. 



where a, b and c are measureable functions on fl with 



> Xy/ac - b 2 I 



1. Introduction 

Elliptic equations in two variables are very well understood, and the regularity 
theory for such equations is significantly stronger than that available for elliptic 
equations in higher dimensions. In particular, Morrey [M] and Nirenberg [N] proved 
Holder estimates for the first derivatives of solutions of uniformly elliptic equations 
in two variables, depending only on bounds for the coefficients: 

Theorem 1. Let fl c M 2 , and set dn(z) = d(z,dCl) for all z € CI. Let u be a 
bounded C 2 (fl) solution of 

(1) aU XX + 2bU X y + CUyy = f 

a b 
b c 

S = \fac — b 2 > everywhere. Then for any a € (0, y/X) there exists C = C(X,a) 
such that for all points p ^ q in fl with d = Tmn{dn(p), dn(q)} > 0, 

\ D <f~ D ^\ < Cd-« sup (|^| + dnf-Vl) ■ 
\p-q\ a a 

These estimates can also be applied to fully nonlinear uniformly elliptic equations 
in two variables, to give Holder estimates for second derivatives, in this case the 
equations have the form 

(2) F[u] = F(D 2 u, Du, u,x)=0 
where u : fl C M 2 -> K, and F : § 2 

x M. 2 x R x 1] — > K is Lipschitz in all variables 
(here §2 is the space of symmetric 2x2 matrices) and uniformly monotone in the 
first argument, so that there exist constants A > A > such that 

(3) XI < [F i] ] < AI 
where F« = 9F {r ' p ' z ' x) . 

or a 
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Theorem 2. Let u e C 3 (fl) satisfy F[u] — in ft C R 2 , and suppose (2) holds. 
Then there exists a — a (A/ A) suc/i that for any points p ^ q in Q with d — 
mm{dn(p),dn(q)} > 0, 

\dMp) - DM<i)\ < Cd - a {]DM m 1} 

where C depends only on A/A and sup |-D-F|. 

In contrast, the situation in higher dimensions is much worse: For fully nonlinear 
equations there is no Holder estimate known for the second derivatives of solutions, 
unless the equation satisfies a concavity condition with respect to the components 
of the second derivatives. The best result available is the following, due to Evans 
[El-2] and Krylov [Kr] (I follow the treatment in [GT]): 

Theorem 3. Let u e C 4 (f2) satisfy F[u] = in f2 where F is a C 2 function of the 
form (2) which is uniformly elliptic (so that XI < F lJ < AI for some A > A > 0) 
and concave with respect to the first argument. Then for any fl' CC 0, 

\D 2 u(p)-D 2 u(q)\ 

blip j j \ o 

P ,qeci' \P - Q\ a 

where a depends on n, A and A, and C depends on n, X, A, |u|c 2 (o) an d d(£l', dfl), 
and on bounds for the first and second derivatives of F (other than the second 
derivative in the first argument). 

In the parabolic case, there are results similar to Theorem 3 (due to Krylov [Kr]): 

Theorem 4. Let u e C 4 (fi x (0,T]) satisfy 

du , 9 . 
— = F(D z u, Du, u, x, t) 
at 

where F is C 2 , XI < [i™] < AI for some < A < A, and F is concave in the first 
argument. Then for any t > and 17' CC fi, 



8 ,t€[rS,,en' V \ P -q\ a + \s-t\»/2 ' |p_g|a + | a _ t |a/2 

\Du(p,t)-Du(p,s)\ 
+ P en'Tl,<T \s-t\(^)/2 ~ ° 

where a depends on n, X and A, and C depends on n, X, A, sup^x^ ^ |£> 2 u|, 
su Pf2x(o.T] \9tu\, d(fl',dfl), t and bounds for the first and second derivatives of F 
(other than the second derivative in the first argument). 

The Morrey and Nircnberg estimates rely either on quasiconformal mapping esti- 
mates or on the fact that in two dimensions the first derivatives of solutions satisfy 
divergence-form elliptic equations. Neither of these methods seems to generalise 
readily to the parabolic setting. There are, however, special estimates known for 
parabolic equations in one space variable, due largely to Kruzhkov [Kz] . 

In this paper I will prove an analogue of Theorem 2 for parabolic equations in 
two space variables. I do not know whether an analogue of Theorem 1 holds. 



\D 2 u(p,t) - D 2 u{q,t)\ , \d t u(p,t) - d t u{q,t)\ 
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Theorem 5. Let f2 be a domain in R 2 . Let u C C 3 (fl x (0, T]) be a solution of the 
fully nonlinear equation 

— = F(D 2 u 1 Du, u, x, t) 
at 

where F is Lipschitz in all arguments and uniformly monotone in the first argument, 
so that XI < [i™] < AI for some A > A > 0. Then for any r C (0, T) and fl' CC n, 

f \D 2 u(p,t)-D 2 u(q,t)\ \d t u(p,t) - d t u(q,t)\ 
., t €[X 6 nA \p-q\« + \s-t\"/ 2 + \p-q\* + \s-t\a/ 2 

\Du(p,t)-Du(p,s)\ 
+ P e^l,t<T \s-t\W ~ ° 

where a depends on X and A, and C depends on X, A, sup nx( - T j (|£> 2 u| + |-Du|), 
su Pfix(o,T] \dtu\, d(Q',dfl), t, and bounds for the first derivatives of F. 

This gives a result of similar strength to Theorem 2 for fully nonlinear parabolic 
equations in two space variables. 

I will also apply similar ideas to parabolic equations in higher dimensions, to 
relax the requirement of concavity to allow just convexity of level sets. 

Theorem 6. Let fl be a domain in M™. Suppose u € C 4 (f2 x (0, T]) satisfies 

— = F(D 2 u, Du, u, x, t) 
at 

where F is C 2 and XI < [i™] < AI for some A > A > 0, and F l ^ kl M l:j M kl < for 
all matrices [M i3 ] for which F tJ M y - = 0. Then for any Q' CC fl and r C (0, T), 

f \D 2 u(p,t) - D 2 u(q,t)\ \d t u(p,t)-d t u(q,t)\ 
s ,te[r 8 Tl P , q en\ \p - q\ a + \s - t\ a / 2 \p - q\<* + \s - t\ a l 2 

\Du(x,t)- Du(x,s)\ 

~ P SUp : -j— r-r- ^ (_/ 

xf=tV,T<.,t<T |s-t|( 1 +«)/ 2 

where a depends on n, X and A, and C depends on X, A, supj 2x ( T ^ (\D 2 u\ + \Du\), 
su Pnx(o,T] \°tu\, d(Q',dfl), t, K, and bounds for first and second derivatives of F. 

In the elliptic case, this extension is trivial because the equation F = is the 
same as the equation tp(F) = if ip is an increasing function with V'(O) = 0. 
However, in the parabolic case there is a big difference between the two equations 

8u 

and 

If F is concave in D 2 u and ip is increasing, then the latter equation is covered by 
Theorem 6 but not in general by Theorem 4. 

I would like to thank Neil Trudinger and Craig Evans for useful suggestions. 
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2. Some background results 



In this section I will recall some results and notation to be used in the proofs of 
Theorems 5 and 6. 

2.1 Function spaces. 

Let f2 be a domain in R™. As usual I will denote by C k ' a the space of functions 
which have all derivatives up to order k Holder-continuous of with exponent a £ 



In the parabolic setting, where we are considering solutions on a space-time 
region Q = Ox (0, T], it is useful to introduce the spaces P k,a {Q) consisting of 
functions u on Q which have the following norm bounded: 



2.2 The Krylov-Safonov Holder estimate. 

The Krylov-Safonov Harnack inequality, first proved in [KS], provides an oscilla- 
tion estimate for solutions of elliptic or parabolic equations. The following Holder 
estimate is a consequence of this. 

Theorem 7. Denote by Q r the region B r x [— r 2 , 0] C R™ x R. Let u : Q R — > R be 
a smooth solution of an equation of the form 



where the coefficients are bounded and measureable, and XI < [a*- 7 '] < AI. Then for 
some C > and a G (0, 1) depending only on n, A/A, and bounds for coefficients, 



2.3 A characterisation of C 1 '" functions. 

I will make use of the following characterisation of C l a functions in terms of 
difference quotients (see [T]): 

Theorem 8. Let u be a smooth function on Br C 1", and suppose there exists a 
constant C such that for all £ <G S 11 ' 1 , h > 0, and x <G B R _h, 



(0,1]. 



U\ P k, a(Q) 




\d b t D a u(x,t)~d b t D a u{y,s)\ 
\x — y\ a + \s — t\% 

\d b D a u(x 1 t) - d b t D a u{x,s)\ 

■Q \s~t\ ^ 



— = a l3 D t D jU + b l D t u + cu + f 



u \p°- a (QR/2) ^ C (Mp°(Qh) + I/U~(Qr)) • 



u(x + h£) + u(x - h4) - 2u(x)\ < Ch 1+ ' 



a 



where a G (0, 1]. Then 



sup 



\Du(x) - Du(y)\ 
\x~y\ a 



< C{a)C . 
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3. TWO SPACE VARIABLES 

Assume (by rescaling if necessary) that u is a C 3 solution on the domain Qi = 
Bi(0) x (—1, 0] C R 2 x R of a fully nonlinear parabolic equation 

^ = F(D 2 u, Du, u, x, t) 

where F is defined on a convex set 5 in S2 x I 2 x I x ^ x (-1,0] containing 
{Ju = (D 2 u(x,t) 7 Du(x,t),u(x,t),x,t) : (x,t) e Qi}- Assume that F is Lipschitz 
on <S>, and satisfies 

AI < [F ij ] < AI 

at each point of S, for some A > A > 0. We will obtain estimates on regions 
Q r = B r (Q) x (-r 2 , 0] for suitably small r. 

3.1 Regularity of the time derivative. 

Let r € (0, 1), and define v T (x,t) — - (u(x,t) — u(x,t — r)). Then v T satisfies a 
parabolic equation onBi x (r — 1,0]: 

dv T _ F(Ju(x, t)) - F(Ju(x, t - r)) 
dt t 
1 f 1 

= - I DF\ sJu{x j) +{1 _ s ) Ju{x ^ T) ■ (Ju(x, t) - Ju(x, t - t)) ds 
T Jo 

= a % i DiDjV T + b l DiV T + cv T + f. 

Here, writing J(s) = sj^x^) + (1 — s)Ju(x 1 t — r), the coefficients are given by 
a ii = £ pi\ J{a) ds , V = FP'\ J(S) ds, c = Jo 1 F'\ J(8) ds and / = %\ J{s) ds. 
Note that AI < o y < AI, and \b l \ + |c| + |/| < CLip(F). 

Theorem 7 applies to give an oscillation estimate for v T independent of t, and 
hence also for u t . In particular, u t is Holder continuous in both space and time, 
and for (xi,ti) and (0:2^2) hi Q1/2 



\u t {x 1 ,ti)-u t (x2,t 2 )\ < C (\x 2 - xxl 2 + \t 2 - h\) a/2 



dF 



dt 



£°°(Qi) 



3.2 Spatial regularity of second space derivatives. 

This is the key estimate: When restricted to each time slice, the function u 
has Holder-continuous second derivatives. This follows by combining the Morrey- 
Nircnberg estimates for elliptic equations with either a perturbation argument or 
an argument using difference quotients. I present the latter argument. 

Fix t e (—1/4,0]. Write w(x) = u(x,t), and (f>(x) = u t (x,t). Then the following 
elliptic equation holds: 

G[D 2 w(x), Dw(x), w(x), x] = 4>{x), 

where G[r, p, z, x] — F[r, p, z, x, t] and <j>(x) — Ut(x,t). Note that AI < G < AI, 
Lip(G) < Lip(F), and <j) satisfies the oscillation estimate derived in Section 3.1. 

Let £ be a unit vector, and for h > let ShW be the difference quotient of w in 
the direction £ with step h: 

5hw{x) = h^ 1 (w(x + h£) — w(x)) . 
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8hW is defined on £>i_/j(0), and satisfies an elliptic equation: Denote Vw(x) = 
(D 2 w(x), Dw(x), w(x), x). Then 

= hr 1 (G[Dw(x + h£)] - G[Vw{x)]) - 8 h <j) 

= h' 1 [ DG\ v{s) ■ (Dw(x + £h) - Vw{x)) ds - 5 h cf> 
Jo 

= a 13 DiDjShW + ip 

Here T>(s) = sVw(x + h£) + (1 — s)Vw(x) and <p — b l Di8hW + cShW + / — 5h<f>, 
where / = f J* 0|p( s ) ds, a" = J* G ij \ v{a) ds, b' = G p% \ v{s) ds, and 2 = 
Jo 1 G z \v (a) ds. It follows that A|w| 2 < a lj ViVj < A\v\ 2 for all v ^ 0, and \<p(x)\ < 
CLip(F)(l + |u| P 2 (Qr) ) + C^-^Lip^) + \u\ P 2 (Ql) ) for x G B 1/i and h < 1/4. 

The Morrey-Nirenberg estimates of Theorem 1 apply to the above equation to 
give the following estimate: If x G B 1 / S (Q) and h < r < 1/16, then 

\DS h w(x + hO - DS h w(x)\ < C(h/r) a (\u\ P 2 {Qr) + r\<p\) . 

where C depends on A/A. In particular, if we choose r = hP for /? G (0, 1) then we 
have for x G £>!/§ and h < 2~ 4 / /3 

|5 fc ff fcJ Dto| < Ch-\l + |w| p2(Ql) ) (V^ 1 ^ + /^-(i-Wi-a)) . 
If we choose 1 > /? > max{0, ^-^7} then we find 

l^h^l <^- 1+£ (l + | U | p2(Ql) ), 

for some e G (0, 1] depending on a and (3, and it follows from the characterisation of 
Qi,a f unc tions in section 2.4 that the second derivatives of w are Holder-continuous 
with exponent s, and 

\w\c*.'(B 1/a ) < C(l + Mp2 (Ql) ) 

where C depends on A and Lip(F). This proves the required spatial regularity of 
spatial second derivatives on the region Q\/%. 

3.3 Time regularity of first space derivatives. 

The spatial C 2,Q estimate established in the previous section can be used to 
deduce an estimate on the continuity of first spatial derivatives in time, using the 
parabolic maximum principle. 

Let £ be a unit vector. The function ShU defined by Shu(x) = {u{x+h£) — u(x))/h 
satisfies a useful evolution equation: 

(3.3.1) -|<5 /lU = a 4 -'A J DA u + ^ 

where supg i/2 \ip\ < C(Lip(F), \u\ P 2(q 1 )). The result of Section 3.2 shows that ShU 
is C l a (uniformly in h), so that 

\S h u{z', t) - S h u(z, t) - DS h u(z, t) .( z '-z)\< C\z' - z\ 1+a 
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on the region Qi/s-h- Young's inequality gives the estimate 

\z'~z\ 1+a < e + Ce'^ \z' -z\ 2 , 
for any e > 0, and therefore 

8 h u(z', t) < S h u(z, t) + D5 h u{z, t) ■ {z! - z) + e + Ce~^ \z' - z\ 2 . 
By the bounds on a y and ip, the function 

t') = S h u{z, t) + DS h u(z, t)-(z' - z) + e + Ce~^ \z - z\ 2 
+ (sup \ip\ + 4CAe-^)(<' - t) 
is a supersolution of Equation (3.3.1) on Qi/i 6 provided h < 1/32, and if z e B 1 / 32 

l+a 

and e < C(A, a)|u|p 2 1 ^ i -j then this supersolution is above Shu on the boundary 
dB x / 16 x ( — 1/256,0]. Therefore by the parabolic maximum principle 5hU is bounded 
by After evaluating at z' — z, and optimizing in e, this becomes 

5 h u{z,t') < 5 h u{z,t) + C(t' - t) 1 ^ + C\u\ P 2 {Ql) \t' - t\. 

A similar estimate from below follows by comparison with a suitable subsolution. 
The desired continuity in time of the first spatial derivatives in Q1/32 follows on 
sending h to zero. 

3.4 Time regularity of second space derivatives. 

The proof of the parabolic C 2,a estimate can now be completed by deducing 
appropriate continuity of the second spatial derivatives in time. We will deduce 
this as a consequence of the previous two estimates. 

On Qi/32 we have the estimates 

\D 2 u(x,t) - D 2 u(y,t)\ <C\x-y\ a 

and 

\Du(x,t) -Du(x,t-r)\ < C\t\^t 

provided t - r > -R 2 /A and x, y e Br/2- Fix x e B 1/6i and s,t e (-1/256,0]. 
Let £ be an arbitrary unit vector. Then we have (provided x + h£ e -B1/32) 

D i D s u(x,t)= D Mx + ^t)-Dsu(x,t) l j\ DiDiu{x+rLt) _ D ^ Diu{x t))dr 

h hj 

and so 

f h 

\D^D^u(x, t) - D 6 D^u(x,s)\ < h' 1 / \(D^D £ u(x + r£, t) - D^D^x, t)\ dr 

Jo 

f h 

+ h' 1 / \(D^D ( u(x + r£, s) - D^D^u{x, s)\ dr 
Jo 

+ h^ 1 \D£u(x + h£, t) - D^u{x + h£, s)\ 

+ h' 1 \D^u(x,t) - D/:u(x, s)\ 

< Ch a + Ch- 1 \s-t\^ 2L . 

Since \s — t\ < 1/256 and x € Bi/64, we can safely choose h = \s — t\ x / 2 < 1/64 and 
still ensure that x + h£ G -61/32- This choice gives 

\D^D ( u(x, t) - D^u(x, s)\ < C\s - t\ a/2 , 
proving the desired estimate. 
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4. Higher dimensions 

In higher dimensions, a similar argument can be used to extend the class of fully 
nonlinear parabolic equations for which second derivative Holder estimates hold: 
Instead of concavity as a function of the second derivatives, it suffices that the level 
sets of F(r,p, z, x, t) as a function of r (for fixed p, z, x and t) are convex. 

Holder regularity of u t follows exactly as in Section 3.1. The key step is to show 
spatial regularity of the second space derivatives in analogy with Section 3.2. 

The argument in Section 3.1 gives the estimate 

\ U t\p°-«(Q 1/2 ) ^ C ( Li P( F ) + Ml~(Q1)) ' 

As in Section 3.2, for any fixed t G (—1/4,0] the uniformly elliptic equation 
G[D 2 w(x) 1 Dw(x), w(x), x] — <p(x) holds, where w(x) — u(x,t), G[r,p,z,x] — 
F[r,p,z,x,t] and <j)(x) = u t (x,t). Note that AI < G < AI and Lip(G) < Lip(F). 
Since the level sets of F in the first argument are convex, and F is uniformly 
monotone and C 2 , there exists a constant K such that 

F klmn M kl M mn < KF kl F mn M kl M mn 

for any symmetric matrix M. But then G = — exp(—KG) is uniformly mono- 
tone and concave in the first argument (with ellipticity constants depending on 
su PQi l^ 2u l)' an d we have G[D 2 w, Dw, w,x] = <f> where = — exp(— K<j>) is Holder- 
continuous on B 1 / 2 - A perturbation result (see [C], Theorem 3) then implies the 
estimate \D 2 w(y) — D 2 w{x)\ < C\y — x\ a for x and y in B\/ & . The arguments of 
Sections 3.3 and 3.4 now apply unchanged to give the full result. 
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